Generation of a superposition of odd photon number states for quantum information 

networks 
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We report on the experimental observation of quantum-network-compatible light described by 
a non-positive Wigner function. The state is generated by photon subtraction from a squeezed 
vacuum state produced by a continuous wave optical parametric amplifier. Ideally, the state is a 
coherent superposition of odd photon number states, closely resembling a superposition of weak 
coherent states (a Schrodinger cat), with the leading contribution from a single photon state in the 
low parametric gain limit. Light is generated with about 10,000 and more events per second in a 
nearly perfect spatial mode with a Fourier-limited frequency bandwidth which matches well atomic 
quantum memory requirements. The generated state of light is an excellent input state for testing 
quantum memories, quantum repeaters and linear optics quantum computers. 

PACS numbers: 03.65.Wj; 03.67.-a; 42.50.Dv 
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Generation of non-classical, non-Gaussian states of 
light compatible with atomic quantum memory has been 
an outstanding challenge driven by various applications 
in quantum information processing. Storage of single 
photon states in atomic memories is at the heart of pro- 
posals for quantum repeaters whereas Fock states 
or Schrodinger cat states, stored off-line, form the basis 
of linear optics quantum computing 0, • Various ap- 
proaches to generation of single photon states compatible 
with atoms have been pursued: single atoms in free space 
£| and in high-finesse cavities 0, |g , and atomic ensem- 
bles 0, IE EJ , and non-classical features such as photon 
antibunching and violation of classical inequalities have 
been demonstrated. However, the dominant contribution 
to the detected state of light in all of those experiments 
has been a vacuum state, due to either high losses or 
high mode impurity. High fidelity single photon states 
[ToL ITU have been reported previously but those were 
not compatible with atomic memories. Schrodinger cat 
states have been demonstrated for microwave radiation 
[iH but not for light. 

In this Letter we report on generation of light, which is 
compatible with atomic quantum memories and which is 
characterized by negative values of the Wigner function. 
The generated states are close to cat states with small co- 
herent amplitude, and converge to a single photon state 
for certain experimental parameters. The compatibility 
with quantum memories is ensured by the purity of the 
spatial, temporal modes of the field. The generated field 
has a well defined Fourier-limited spectrum with a band- 
width around 9 MHz in the low gain regime, is frequency 
tunable, and has a perfect spatial Gaussian profile. 

We follow the essence of the proposal 0] m sub- 
tracting a photon from squeezed vacuum by transmit- 
ting squeezed light through a weakly reflecting beam 
splitter and detecting a photon in the reflected beam. 
In 0] this idea was used for demonstration of gener- 
ation of non-Gaussian states. In the absence of losses, 
the state in case of, e.g., 3 dB of squeezed input is 



* ~ |1) + 0.29|3]_+- 0.09J5) + ... in the Fock basis. It has 
been shown in 0, ITU ITr3 | that such a state well approx- 
imates a small-amplitude Schrodinger cat state. 

In order to generate a narrow-band state of light in a 
pure spatial mode, we make use of a continuous wave (cw) 
optical parametric oscillator (OPO) below threshold for 
generation of the initial squeezed field. Such a device has 
been successfully used for generation of squeezed and en- 
tangled Gaussian states compatible with atomic targets 
[H TE1 111 HI HT|. In an OPO the non-linear optical 
medium is placed inside an optical resonator (Fig. |T]l. 
thus providing light generated within well defined fre- 
quency and spatial modes. At the same time the res- 
onator significantly enhances the spectral brightness of 
the selected modes. All those features prove advanta- 
geous for the OPO-based source compared to single-pass 
sources 0, 0, ^| as far as compatibility with atomic 
memories is concerned. The cw character of the exper- 
iment requires a multimode theory for its adequate de- 
scription. We have developed this theory based on the 
formalism of Gaussian states and Gaussian covariance 
matrices [53] (see also an alternative theoretical multi- 
mode analysis for a different detection model |23|). 

In the experiment (Fig. ^| the output of the 
TkSapphire laser around 852 nm (close to the Cesium 
D2 transition) is frequency doubled to generate up to 
150 mW of the blue OPO pump beam (426 nm), while 
a small part of the laser beam is split off to act as the 
local oscillator (LO). The OPO is an 81 cm long bow-tie 
cavity with a 10 mm, temperature-tuned KNb03 crys- 
tal. The down-conversion is Type I with 90° non-critical 
phase matching, hence all down-converted photon pairs 
are generated in the same polarization and in the same 
direction defined by the spatial mode of the resonator. 
The OPO cavity is locked on resonance with the laser 
frequency, so that the degenerate down-conversion mode 
is efficiently generated. The output coupling mirror of 
the OPO has a transmission of 12.7% and the intracav- 
ity losses are 1.0%. In addition to this, the blue pump 



2 




HD 



T 



Scope 



i 



APD IF 



ffl A 










G 


= 2.3 




CO 




n, 


= 0.56 


'O 2 
c 

15 








i_ 

-t-j 


i-requency [Mi-izr~-*-*— 


u 

a 

W -2 


2 4 



FIG. 1: (Color online) Experimental setup. SHG - second 
harmonic generator, OI - optical isolator (to prevent scattered 
LO light from entering the OPO), A/2 - half wave plate, PBS 
- "magic" beam splitter, PZT - piezo mounted mirror, HD - 
homodyne detector, IF - interference filter, APD - avalanche 
photo diode. The inset shows the squeezing spectrum ob- 
tained from a Fourier analysis of the raw data. The blue 
curve is the best fit of the theoretical spectrum. The fitted 
parameters (OPO gain 2.3, and total efficiency, rjt = 0.56) 
match well the experimentally estimated values. 



induces around 1% absorption of the IR cavity field due 
to the light-induced loss process in the crystal - BLIIRA 
[24|. Altogether, this gives a cavity bandwidth of 4.4 
MHz HWHM and the output coupling efficiency of the 
cavity around w npn = 0.86. More details of this setup 
can be found in [H H3 . 

In order to characterize squeezing generated by the 
OPO and the losses in the homodyne channel, we mea- 
sure the squeezing spectrum of the OPO. The OPO out- 
put is mixed on a balanced 50/50 beamsplitter with a 
few-milliwatt LO beam. The differential photocurrent is 
measured by two low noise photodetectors and its Fourier 
transform is digitally obtained. The results are presented 
in the inset in Fig. Note that squeezing is observed 
across the entire theoretical bandwidth of the OPO (with 
the exception of frequencies below 30 kHz), and not just 
within a narrow frequency band. This is necessary for 
the photon subtracted state generation. From this data 
we infer the OPO gain and the overall losses. The overall 
propagation efficiency 0.66 coincides reasonably well with 
V = VoPoVpriVhom) 2 calculated from the measured val- 
ues for rjopo, propagation efficiency from the OPO to de- 
tectors ri pr — 0.88 and homodyne efficiency r] hom = 0.96. 
Our homodyne detectors have the quantum efficiency of 
0.90 which along with the small dark current lead to the 
effective detector efficiency rjd = 0.85, giving the total ef- 
ficiency of detection rjt = 0.60 (close to r\t — 0.56 inferred 
from squeezing data). 

For the photon subtraction experiment the down- 
converted output of the OPO is directed onto a "magic" 



beam splitter (polarizing beam splitter cube in combina- 
tion with a half-wave plate) which reflects a small por- 
tion (3-7%) of the light towards a photon counter - APD 
(Perkin-Elmer SPCM-AQR-13, dark count rate 160 s" 1 ). 
Note that for low reflectivity of the "magic" beam splitter 
the probability of a two-photon click is very low, i.e., a 
click of the APD almost for sure corresponds to the sin- 
gle photon detection. Besides the degenerate squeezed 
mode the OPO generates thousands of, non-degenerate, 
modes which are also phase matched and resonant with 
the cavity. Photons in those modes produce unwanted 
counts in the photon subtraction channel, and therefore 
have to be filtered out. Towards this end, on the way to 
the APD, the light is filtered with three consecutive cav- 
ities and a spectral filter which block all the OPO output 
modes except for the degenerate squeezed mode (Fig. ^| . 
The OPO cavity and the three filter cavities are tuned 
in resonance with the degenerate OPO mode by inject- 
ing a weak coherent seed beam into the OPO. The count 
rate of the APD is then used to produce the error signal 
for the frequency locking. A Lab View program monitors 
this count rate while it scans all cavities at different low 
frequency rates until the program finds a setting for the 
four piezo control voltages which gives an expected max- 
imal count rate. Under stable lab conditions the lock 
is reliable down to a count rate of approximately 5000 
s -1 . With the cavities lock points established, we in- 
ject the blue pump into the OPO and turn down the co- 
herent seed beam keeping the lock operational now with 
only the down-converted light. By turning the "magic" 
beamsplitter to full reflection, we measure the maximum 
count rate as a reference for the subsequent setting of the 
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FIG. 2: Part of a typical trace obtained after application of 
the temporal mode function (shown in the inset) to each of 
the 20,000 raw data segments. Each segment of raw data is 
multiplied by the mode function and integrated, producing 
one phase/quadrature point for the trace. The phase of the 
LO has been scanned at roughly l7r/s. The temporal mode 
function is shifted with respect to the APD click due to the 
delay in electronics. 
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FIG. 3: (Color online) Figures a,b) show the density matrices and Wigner functions (including a top-down contour plot) 
reconstructed from experimental data, a) the OPO gain 1.8 b) gain 2.3. The presented plots are averages of 50 separate states, 
each reconstructed from a 20,000 points quadrature trace as explained in the text and corrected for the 85% detector efficiency, 
c) state calculated from full multi-mode theory with the same temporal mode function and the same overall propagation 
efficiency r\ = 0.64 as for the experimental state b). d) state calculated with perfect efficiency r) = 1. 



beamsplitter at the chosen reflection ratio. At the same 
time we get an estimate for the total production rate of 
photons into the degenerate mode (around 2* 10 6 s _1 for 
the parametric gain 2.3). 

The LO phase is scanned to perform the tomography 
of the state. For every click of the APD, a data segment 
of 1000 points is recorded within a period starting 1 /is 
before and ending 1 /_ts after the trigger event. 20,000 
such segments are stored and comprise one full measure- 
ment of the different quadratures of the quantum state. 
The total acquisition time of a full measurement is only 
around 3 seconds and a new measurement can be initi- 
ated a few seconds later. The high spectral brightness 
of our source thus allows faster data acquisition than in 
single pass experiments where it may take hours. 

During the initial data processing each 2 /is segment 
of experimental data is turned into a single quadrature 
measurement using a particular temporal mode function. 
If the OPO is in the intermediate to high gain regime, 
the field populates many modes, and the optimum mode 
function should be chosen as the one populated by a 
quantum state with the strongest non-classical features, 
e.g., the one with the largest negative value attained by 
the Wigner function. A fairly simple Ansatz for the mode 
function k e — 27r 7l<l — ^ e -2*K.\t\ ^ which j s the Fourier trans- 
form of the product of two Lorentzians with half-widths 
2-7T7 and 2nn, shows close to the optimal performance. 
The optimal 7 and n are close to 2ir * 9 MHz (OPO 
width) and 2tt * 48 MHz (filter cavity width). The tem- 
poral shape of the generated state is shown in Fig. 
Every 2 /js/1000 points data segment is then multiplied 
by the mode function. This gives a number which is the 
measured quadrature amplitude value for that segment. 
This procedure is equivalent to the use of a pulsed lo- 
cal oscillator with the temporal shape of the mode func- 



tion. A full measurement set consists of 20,000 quadra- 
ture values, along with their corresponding phase values 
(obtained from the time of their acquisition via an inde- 
pendent phase calibration). A data set obtained in this 
way is presented in Fig. |2J The lack of points around 
zero values of the quadrature phase amplitudes is a clear 
signature of the non-Gaussian character of the radiation. 
We perform the same procedure for the vacuum state, 
and use the result for normalization. 

To reconstruct the density matrix from the data we use 
the iterative maximum- likelihood method of |25| . based 
on [2^| - for reviews of this and other methods, see |27| . 
In order to obtain the Wigner function of the optical 
field irrespectively of the detector inefficiency we cor- 
rect the density matrix for the detector efficiency us- 
ing the algorithm of [25j. The density matrices and the 
Wigner functions of the photon subtracted squeezed vac- 
uum states for two different degrees of squeezing are 
shown in Fig. [3^,b. Each presented state is an aver- 
age of 50 individually acquired 20,000-point measure- 
ments. The averaged values of the Wigner functions at 
the origin are, respectively, W(0, 0) = -0.028 ± 0.010 
and W(0, 0) = -0.040 ± 0.014 (we note that even the 
states reconstructed without correction for detector effi- 
ciency have small negative dips) . The negative regions of 
the Wigner functions show their profoundly non-classical 
character. The Wigner functions are elongated along 
the direction of the parent squeezed states signifying the 
presence of three photon and five photon contributions. 
In the limit of very low degree of squeezing the elongated 
shape would be reduced to a round shape characteristic 
for the single photon state. Fig. [3]; shows the state calcu- 
lated following our theoretical model with the efficiency 
7/ = 0.64 corresponding to the experimental efficiency of 
the state shown in Fig. ^jp. This number is a product of 
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•q = 0.66 obtained from the squeezing data and an extra 
factor 0.97 which partially accounts for dark counts of 
the APD. As seen from Fig. 13b and Oh, theory and ex- 
periment are in good agreement. The small discrepancy 
may be attributed to an admixture of other modes in the 
trigger channel. 

Due to the cw character of the squeezing produced in 
the experiment, the single mode theory of |13lll5lll6ll28| 
is not fully applicable. We have developed an exact the- 
oretical treatment, which is possible because the output 
from the OPO is a Gaussian state fully described by its 
second order moments, i.e., by the two-time correlation 
functions of the output (&t (*)&(*')) and (b(t)b(t')) 0. If 
the transmission (reflection) of the "magic" beam split- 
ter is |t| 2 (\p\ 2 = 1 — | -r | 2 ) and if the filters for the trig- 
ger field can be modeled by a single Lorentzian of half 
width 2ttk, the annihilation operator for the mode, caus- 
ing the click of the APD at t = t c can be expressed 
as a te = 2p v ^7TK J_ e oo exp(-27rK(t c - t))b(t)dt + F noisc 

where -Fnoisc represents vacuum field contributions from 
the "magic" beam splitter and from the frequency filters, 
and its precise form does not need to be specified. For 
an arbitrary mode function u(t), we obtain a similar ex- 
pression for the annihilation operator of the field mode 
extracted by the corresponding LO detection within the 
time mode u(t), a u — t f_ u(t)b(t)dt + G no i SC , where 

G n oisc is a vacuum noise term. Losses and detector effi- 
ciencies are incorporated by reduction factors on the first 
terms and additional vacuum noise contributions in the 
equations for &t c and a u . The second moments of the 
four field quadrature variables for these operators are di- 
rectly obtained from the output correlation functions and 
the filter and mode functions. Defining the 4x4 covari- 
ance matrix 7 by its elements 7^ = 2Rc(yiyj) , with the 
four Hcrmitian field quadratures forming a column vector 
y = (yi, 7/2, 2/3, 2/4) T , the Wigner function can be written 
explicitly, W(y) = l/(V d et(7)) exp(-y T 7 - 1 y)- This 
Gaussian Wigner function fully characterizes the quan- 
tum state of the trigger and LO modes prior to the APD 
detection event. The action of a photo detection event 
is described by application of the corresponding field an- 
nihilation operator to the quantum state, which in the 



Wigner function representation is represented by a com- 
bination of differentiation and multiplication by the yi 
arguments 29]. The reduced state associated with the 
homodyned field mode is subsequently obtained by trac- 
ing over the trigger field state, i.e., by integrating over 
the trigger mode quadrature arguments in W. These op- 
erations can be done explicitly, and the resulting single 
mode phase space function is readily obtained (Fig. [3k). 
For further discussions of this treatment, see |22j |. 

The state generated by photon subtraction from single 
mode squeezed vacuum is a superposition of odd- number 
Fock states 0, which overlaps with a cat state \a) — \— a) 
with better than 99% fidelity for coherent amplitudes less 
than a — 1.2. In the present cw experiment we have sev- 
eral modes populated, hence we are bound to get contri- 
butions from even number states, and e.g. the vacuum, 
which grow with the degree of squeezing. The states 
shown in Fig. |2K,b overlap with cat states with coher- 
ent amplitudes a = 1.05 and a = 1.3 with a fidelity of 
53%. For our squeezing, the overlap between the lossless 
state (Fig. |3t) and a cat state would be still quite good, 
0.95 and 0.90, respectively. A major improvement of the 
purity of cat states can be achieved by reducing optical 
losses. It seems feasible to obtain an OPO escape effi- 
ciency rjopo — 0.95, a propagation efficiency rj pr — 0.97 
and a homodyne efficiency rjhom — 0.99, thus reaching a 
value close to 0.87 for the total efficiency. Even higher 
purity can be obtained with purification procedures pro- 
posed in 0,0. 

In conclusion we have demonstrated experimentally a 
new type of non-classical state of light, the photon sub- 
tracted squeezed vacuum state based on an optical para- 
metric oscillator below threshold. The results are in good 
agreement with a complete multi-mode theoretical anal- 
ysis. With the time mode function of 100 nsec, and the 
transform limited bandwidth around 10 MHz, frequency 
tunable around Cesium D2 transition, the source is com- 
patible with atomic memories and networks. 

After submission of this paper, generation of an optical 
"kitten" state was reported by Ourjoumtsev et al. |30j . 

This research has been funded by EU grants COV- 
AQIAL and QAP. 
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